Heavier Higgs Particles: Indications from Minimal Supersymmetry 
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We use the most recent data on the Higgs-like resonance h observed at 125 GeV to derive 
information about the mass of the heavier Higgs particles predicted by Minimal Supersymmetry. 
We treat as independent parameters the couplings of h to top quark, beauty and massive vector 
bosons and, in this three dimensional space, we locate the point realizing the best fit to data 
and compare it to the position of the Standard Model point and to the region of coupling values 
accommodating heavier Higgs particles in Minimal Supersymmetry. We conclude that mass values 
320 < Mh 360 GeV are compatible at 2a with the best fit of couplings to present data, larger 
values being compatible at the la level. Values of 1 < tan /3 < 6 are compatible with data. 
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In a recent paper [1| , we confronted Minimal Supersymmetry for the Higgs sector with the particle discovered 

by ATLAS and CMS at 125 GeV. At that time, only the value of the mass was available, which, as noted also 
by other authors [1, H| , is remarkably consistent with Minimal Supersymmetric Standard Model (MSSM) , albeit 
with a value of the scalar top mass around 4 TeV. The presence of a tenuous bump at about 320 GeV in the 
CMS data for the ZZ channel could have been taken as a hint of the heavier, MSSM, + Higgs particle, H . 
With the higher statistics now available there is no trace left of such a structure. 

In this note, we confront MSSM with the most recent data of the resonance /i(125) provided by ATLAS Q, 
CMS 5) and Tevatron (6j. Mass is given as: 

D. A26.0 ± 0.4 ± 0.4, ATLAS\ , , 

rC; Mh - \125.3 ±0.4 ±0.5, CMS ) <- > 

■ For simplicity, we shall use in the computations of cross sections and decay rates the conventional value Mh = 125 
> GeV. 

^sO . The experimental signals for decays are given by the ratios: 

_ fa x BRj) cxpt 

: Ml ~ fa x BRj)sM 1 ' 

ON 

for the channels i = WW, ZZ, bb, rf, 77, for ATLAS Q and CMS @, and i = bb for Tevatron @. 

We aim to obtain limits on the mass of H and of its closeby partners A and if ± and on the tan /3 parameter. 
We find that, to one standard deviation, MSSM is consistent with present observation for Mh ^ 360 GeV, 

■ and for tan (3 > 1. 

We study the problem in the (q, cj,, cw) space of couplings of the h resonance to the top and beauty quarks 
and to the vector bosons, Cj representing the ratios of the MSSM to the Standard Model (SM) couplings. We 
set cz = cw and use the same coupling ct for all up fermions, e.g. charm, and q, for all down fermions, e.g. r. 

Two different types of cross sections for the h production appear into the definition of fii given above, namely 
the gluon fusion cross section cr g f and the vector boson fusion one, <r VBF . These are parametrized according to 

<7 g f = (<7gf) SM cf (3) 
°Vbf = (°Vbf)sm c W (4) 

Cross section ((31) is dominated by the top quark loop and it applies to the first four decay channels listed above. 
The cross section in ((4]) appears in the 77 channel due to the stringent cuts applied by ATLAS and CMS to 
the diphoton distribution so as to enhance the purity of the 77 signal. In this case we write: 

^ = faf) S M c t + ( 5 'vbf)sm c \V ( 5 ) 

where the bar represents the effect of the experimental cuts. The ratio a/asM depends only on the ratio: 

R = ^ (6) 
cr gf 
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In Ref. [l| one finds a detailed calculation of R at 7 TeV using ALP GEN Q supplemented by the cuts specified 
by CMS in Ref. @|. No difference using ATLAS cuts is found. We found R = 2.97 @ 7 TeV, whereas we 
find, with the same tools, R = 3.1 @ 8 TeV. Neglecting such a small deviation, we can use the cross sections 
computed in [lj to compare with the experimental values for /i,-, which are obtained by combining data at 7 
and 8 TeV. 

Widths in Vector- Vector and fermion antifermion decays of h are parametrized by rescaling the Standard 
Model values by the appropriate factors, cfy, etc. As for the 77 width, we write: 



\g icw + gi/2 Ct\ 2 

\gi +gi/2\ 



nrr) - r(77) SM ' , 7r7 ^ 



where 51,1/2 are the loop contributions of W and top [lOj, with: 

5 i2 = 8.32, ff i /2 = -1.84 (8) 

The total width is obtained from the sum of the individual widths, and branching ratios are obtained accord- 
ingly. 

We discuss in the end the possible contribution of scalar top quarks in the gg and 77 loops. We anticipate 
however that the eventual effect of scalar top quarks in the loop simply amounts to a redefinition of Cf , leaving 
the structure of the above equations unchanged [20j |. 

To find the best fit values of the c's, we minimize the usual x 2 function: 

2 _ [/Xq - Gg(Ct,Cb,Cw)} 2 , Q s 

a a 

where fi a runs over the values of /x, Eq. @ given by ATLAS, CMS and Tevatron, a a the corresponding error and 
G a the theoretical expressions introduced above in terms of the couplings. We have in all eleven experimental 
data and three parameters, namely eight degrees of freedom. 

As a first step we identify the points in this space which minimize the y 2 function. Setting conventionally 
c-w ^ 0, we obtain the two solutions shown in Fig[TJ corresponding to c t < and q > 0. The two solutions 
have the same value of x 2 , with x 2 /d.o.f = 0.6. 




FIG. 1: la (inner) and 2a (outer) regions in the space (ct, Ct,, cw) around the best fit points obtained by minimizing 
the x 2 function. The fit includes ATLAS, CMS and Tevatron data, see text. Solutions are at (ct, C&, cw) + = (0- 886, 1.14, 
1.19) and (c t , Cb, cw)~~{-0. 908, 0.642, 0.666). We also add a point for the Standard Model value (green point). 



The relative sign of c t and cw is relevant to the 77 decay, see Eq. ([7]). It is no surprise that there is also a 
best fit solution with ct < 0, particularly in a situation where the central 77 signal exceeds the SM value 0,01 ■ 
A similar situation was also found in recent two-parameter fits 0, Illl - fl3j . 
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The inner and outer shells in Fig. Q] represent the la and 2a regions. In the same figure we place the Standard 
Model point (1, 1, 1), which falls within la from the best fit to data. We will focus on the positive solution. 

Next we study the surface S : (c t (tan/3, m# ), Cb(tan/3, ran ), cw(tan /3, m# )) predicted by Minimal Supersym- 
metry. Here we briefly recall its definition. 

To start, let us recall that in the basis (Hd,H u ), the mass matrix of the CP-even Higgs fields is given by: 

Ml = M%( C ° s2 / R - C0S ^)+Ml( si * 2 / R " cos ^f + °.) (10) 



- cos /3 sin p sin ft J V ~ cos P sm P cos P J \0 $ 

with S the radiative correction 



i = J^ G ,(M,)«l„ g ffi (11) 
7r 2 sm p \ M t I 

The first term in M 2 S arises from the so-called Fayet-Iliopoulos term determined by the gauge interaction. In 
the radiative corrections we have kept only the top-stop contribution. We determine the system of eigenvalues 
and eigenvectors of this matrix expressing Ma in terms of the CP-even Higgs boson mass, Mjj, and the known 
values for Mz and Mh = 125 GeV, according to: 

M 2 A = Ml + M 2 h - M 2 Z - 5 (12) 

The latter relation is obtained by equating the trace of M. 2 S to Mjj + Mi (the sum of eigenvalues). 
In the calculation of eigenvalues and eigenvectors, we also express the quantity t, defined by: 



/ J Mi Mi \ 

'^("t) (13) 

as a function of tan/3 and Mh, i.e. we find t — ^(tan/?, Mjj ). This can be obtained from the eigenvalue 
difference squared, (M H — Ml) 2 . The function F arises as the solution of a second order equation. However, 
the other root is spurious, as it does not reproduce the degeneracy of the heavier Higgs boson masses as 
Mh — > oo, see below, Eqs. (ITBl . More on the function F can be found in 

In the literature, it is customary to use Ma as an independent parameter. Here, Ma is obtained from eq. 
([T2"T) and from the function ^(tan /?, Mh)- In addition: 

M 2 H± = M 2 A + M 2 , (14) 

For future reference, we record some values of Ma and M H ± (in GeV): 

M A = [315,358]; M H ± = [325,367] (M H = 360) 
M A = [600,619]; M H ± = [604,624] (M H = 620) 

tan/3 = [1,6] (15) 
We next proceed to the computation of the normalized eigenvector corresponding to /i(125), namely: 

and define the surface S parameterized by tan/3 and Mh in the (ct,Cb,c\v) space: 



c t (tan/3, M H ) = ^ l + t& f P ^(tan/3, M H ) (17) 
tan p 

c 6 (tan^, M H ) = yi + tan^/3 5hd(tan/3, M H ) (18) 

1 tan B 

c w {tanB,M H ) - . S hd (tfmP,M H ) + . 5 hu (tan/3, M H ) (19) 

VI + tan 2 P V 1 + tan P 

In Fig. [2] we show the la and 2a regions around the best fit solution with ct > 0, intercepted by the surface 
S described above, in the parameter range: 

320 GeV < M H < 620 GeV; 1 < tan /3 < 6 (20) 
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FIG. 2: The intersection of the surface S in the (ct,Cb,cw) space, as predicted by MSSM, with the la and 2a regions 
around the best fit point to ATLAS and CMS data. 



The so called decoupling limit is manifested here in the fact that S heads straight into the Standard Model 
point as the mass Mjj is increased. To be definite, we restrict to the mass range (J20J) - 

A better display of the situation is obtained by fixing cw = 1 and considering the intersections of the ler, 2a 
ellipsoids with that plane. This is illustrated in Fig. [3] where we also plot the surface S as embedded in the 
(ct,Cb) space. 

The (blue) points on the left side of S correspond to Mh — 320, 360 and 620 GeV, the intermediate value 
marking the border of the la region. The (red) points crossing the surface from the left to the right are 
associated to different values of tan/3, from 1, at the lefthand side, to 6, at the righthand one. The best fit 
point projects inside the MSSM surface, but is above it (&w = 1.19) in the three dimensional space. 




c, 



FIG. 3: The surface S as embedded in the (ct, Cf,) space at cw = 1. The la and 2a regions are the intersections of the 
ellipsoids at c t > 0, Figs. 1 and 2 with the plane cw = 1. The outer side and the short inner side of S correspond to 
Mh = 320 and Mh = 620 GeV, respectively. The intersection with the la border corresponds to Mh = 360 GeV. Red 
points mark values of tan/3 from 1 (lefthand side) to 6 (righthand side). The best fit point projects inside the MSSM 
surface, but is above it (cw = 1.19) in the three dimensional space. 

The outcome of this analysis is that, even though we lost the footprints of an heavier Higgs at 320 GeV, 
the hypothesis that the second + particle predicted by Minimal Supersymmetry might just be lurking in the 
background is still viable. 
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To la, the region Mh ^ 360 GeV is all allowed, along with 1 < tan/3 < 6, consistently with the limits 
arising from the absence of Flavor Changing Neutral Current B s decays. Data are compatible with a regime of 
incomplete decoupling, the 'intermediate - coupling regime' of Ref. [3j. However, MSSM could well be realized 
in the full decoupling mode. Consistency to la with SM is also found in [l4| . 

In Ref. [l[ we had reported the values of a x BR for a hypothetical H at Mh = 320 and 500 GeV. We 
close here by giving the picture for Mh at 360 and 620 GeV and comment on production and decay of the 
corresponding pseudoscalar and charged Higgs particles, A and H^. 



b 0.2 
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FIG. 4: Production of H and A at the LHC (total energy 8 TeV) from gluon fusion. Cross sections (pb) are given for 
tan/3 = 1.5, as functions of M^(TeV), in the 2a allowed range. 

Cross sections for H and A production by gluon fusion are shown in Fig. 2] in the full 2a allowed mass range, 
360 < Mh < 620 GeV. Cross sections have been computed with ALPGEN [7j, and corrected with the gluon 
K-factor taken from [To] ]. 

Branching ratios for the dominant decay channels of H are reported in Fig. [51 for Mh = 360 GeV, full lines, 
and Mh — 620 GeV, dashed lines, for 1 < tan/3 < 6. Decays are dominated by the ti channel, with vector 
boson and bb channels still visible in the lower part of the mass range. 

A and 77 ± decay branching ratios can be obtained with HDECAY [16]. In the la allowed mass range, top 
quark channels, ti and tb respectively, are are by far the dominant ones. Production at the LHC and decays of 




FIG. 5: Branching ratios for the dominant deecay channels of H. Solid lines for Mh = 360 GeV, dashed ones for 
Mh = 620 GeV. Decays are dominated by the ti channel, whereas vector boson and bb channels are still visible in the 
lower part of the mass range. 



the H have been extensively studied in Refs. 17 1. The dominant mechanism, in the mass range considered 
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here, is by 6-gluon fusion: 



b + g^t + H- (21) 

followed by H~ -> t + b. 

Finally, we comment on possible effects of scalar top quarks in the hgg and /177 loops, in view of the large 
contribution of scalar tops to the Higgs mass correction, Eq. (fTTj) . Such effects have been considered in [l8|, • 
Unlike the mass correction, the 77 and gg loops are ultraviolet convergent, therefore heavy scalar particles 
decouple in the infinite mass limit. The only interesting case is the contribution of a light stop, the other one 
being correspondingly massive in order to satisfy Eq. (fTTj) . This case has been considered in [l8| and leads, 
neglecting mixing amongs scalar quarks, to: 

A t+t _ t , ( m t\ A (t) _ i 1 t n oc ( m t 



a, 1 + {wJ ;Or 1+0 ' 25 UtJ = MM <' } (22) 

where t = 4M? /m| and r t = 4mf/m|, and A ^/ 2 are defined in [l8|, Gil- As noted before the additional 
contribution simply amounts to a redefinition of Ct and therefore does not affect the best fit procedure. In the 
no-mixing approximation, the tan j3 dependency of the top and scalar top loop amplitudes is the same so that 
the effect of scalar top simply leads to an effective loop t-quark coupling: 

c* oop (tan/3, M H ) = f t {M h ) c t (tan/3, M H ) (23) 

This amounts to a 3% correction for Mj = 500 GeV. 

In conclusion, present experimental information is consistent with the heavier scalar particles of MSSM to 
be in a relatively low energy region, between 320 and 620 GeV. A meaningful search can be done at the LHC, 
even at present luminosity. We underscore that higher Higgs particles still represent one of the most interesting 
test-tables of Minimal Supersymmetry. 
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